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ABSTRACT 

We develop chiral perturbation theory for baryons in quenched QCD. Quenching (the elimina- 
tion of diagrams containing virtual quark loops) is achieved by extending the Lagrangian method of 
Bernard and Golterman, and is implemented in a theory where baryons are treated as fixed velocity 
sources. Our method requires that the octet baryons be represented by a three index tensor rather 
than by the usual matrix field. We calculate the leading non-analytic corrections to the masses 
of octet and decuplet baryons. In QCD these are proportional to M%. We find that quenching 
alters the terms, but does not completely remove them. In addition, we find non-analytic 
contributions to baryon masses proportional to M w and M%logM n . These terms, which are arti- 
facts of quenching, dominate over the M% terms for sufficiently small quark masses. This pattern 
of corrections is different from that in most mesonic quantities, where the leading non-analytic 
terms in QCD (proportional to M^logM,,-) are removed by quenching. We also point out various 
pecularities of the quenched theory, most notably that the A baryon can decay (if kinematically 
allowed), in the sense that its two point function will be dominated at long Euclidean times by a 
nucleon plus pion intermediate state. 



1 Introduction and Summary 



A central goal of numerical simulations of lattice QCD is to calculate the hadron spectrum 
and compare with experiment. Agreement between simulations and experiment would pro- 
vide a crucial test that QCD is the correct theory of the strong interactions in the low energy 
domain. Present simulations have yet to achieve this goalfj largely because of various ap- 
proximations that must be made to make the calculations computationally tractable. Most 
important of these is the use of the quenched approximation, in which the fermion determi- 
nant is approximated by a constant, so that there are no internal quark loops. The other 
significant approximations are the use of quarks with masses heavier than their physical 
counterparts, and the use of finite lattice spacing. 

This paper is the first of two in which we study the importance of the quenched approxi- 
mation for baryon masses. In this paper we develop chiral perturbation theory for baryons in 
the quenched approximation and calculate the dominant contributions from one-loop graphs 
to both octet and decuplet baryon masses. In the companion paper we work out some impli- 
cations of these results |§ . In particular we use them to estimate the size of quenching errors 
in baryons, and show how they lead to a better understanding of how to do chiral extrap- 
olations in quenched QCD (QQCD). This knowledge has already been useful for analyzing 
results of simulations M. 

To provide a context for the quenched results, it is useful to recall the general form of 
the chiral expansion of baryon masses in QCD (see for example Refs. j|, |5], |5|), 

M bary = M + c 2 Ml + c 3 M* + c iL M* log M n + c A M^ + . . . . (1) 

Here Mo is the mass in the chiral limit, and M n is a generic pseudo-Goldstone boson mass. 
The dominant contributions come, in fact, from the K and rj mesons. The q are combinations 
of coefficients appearing in the chiral Lagrangian. The expansion consists of terms which are 
analytic in the quark masses — those proportional to M^ oc (m u + m^) and — and non- 
analytic terms proportional to M^ and logM^. The latter arise from infra-red divergences 
in loop graphs, and thus are multiplied by constants (03 and c^l) which are determined in 
terms of leading order chiral coefficients. For example, C3 depends on the pion-nucleon 
couplings F and D (defined precisely below). Thus, if one knows the lowest order chiral 
coefficients, one can predict the form of the leading non-analytic terms. The same is true 
for mesonic quantities (e.g. M^ and f n ), but with baryons the leading non-analytic terms 
contain one less power of M v than the first analytic corrections (M% versus M%), and are thus 
enhanced in the chiral limit, while with mesons (as with the c^l term above) the enhancement 
is only logarithmic. 

In the quenched approximation we find that the expansion is different, 

K*y = M o + * 4M n + 5 c\ L Ml log M n + c\Ml + c\Ml + c\ L M* log M w + c\M* + .... (2) 

The constants are now different combinations of the coefficients in the quenched chiral La- 
grangian. We have calculated the constants c\, c\ Ll cl and c|, all of which are given in terms 
of lowest order coefficients. We have not calculated c\ L , but for sufficiently small M^ (which, 

1 Opinions differ as to how close we are to reaching this goal. For example, Ref. |Q is optimistic while 
Ref. H is more pessimistic. 
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in numerical simulations, is adjustable) the M% terms should be a small correction compared 
to the lower order terms. 

The dimensionless constant 5, which multiplies the c\ and c\ L in Eq. |2], is a pure quenched 
artifact. It appears also in mesonic quantities, and its likely magnitude is S ~ 0.1, although 
there is no definitive determination J7|. In Eq. it multiplies terms which will dominate 
the corrections for small enough M n . If such terms are significant, then we should not 
trust the quenched results. Fortunately, as discussed in the companion paper, the numerical 
evidence suggests that these terms are small, and so are largely curiosities. They do, however, 
give another example of the sickness of the quenched approximation in the chiral limit, a 
phenomenon first pointed out for mesonic quantities in Refs. || ||. 

Of more practical interest is the fact that the non-analytic terms proportional to M% 
survive quenching, albeit multiplied by different constants than in QCD. Naively, one might 
expect meson loops to require the presence of underlying quark loops, and thus that quench- 
ing would remove all the M% terms. This is indeed true for mesonic quantities, where only 
artifacts proportional to 5 survive. Why this is not true for baryonic quantities is explained 
in Sec. |3|. 

The presence of the M% terms has practical implications. When doing chiral extrapo- 
lations to quenched results one should fit to a form including these cubic terms. It turns 
out that certain linear combinations of baryon masses have expansions beginning at O(M^), 
and these combinations can be used, in principle, to extract the constants cf. With these in 
hand, we can then make estimates of the errors in baryon masses due to quenching. The idea 
is to compare the size of the contributions of pion loops in the two theories, and assume that 
any difference is a quenching error. The uncertainty in these error estimates can be reduced 
by forming suitable combinations of baryon masses. Results from a pilot study along these 
lines are presented in the companion paper. 

In present simulations, the splitting between the octet and decuplet baryons is substan- 
tially smaller than the pion mass, Ma — < M n . Thus in our calculations we treat 
Ma — M^ as a small parameter, and expand about the limit that the octet and decuplet 
baryons are degenerate. Eventually, simulations will be done in the opposite limit, in which 
case the A should be treated as heavy and integrated out of the effective theory. As ex- 
plained below, it is straightforward to adapt our results to this new limit or, indeed, to 
any intermediate value of (Ma — Miv)/M 7r . An interesting phenomenon which occurs when 
M n < Ma — Mn is that the quenched A decays. In Sec. we explain why this happens, 
and why it might not have been expected. 

Finally, we note that the methods presented here can be extended straightforwardly to 
"partially quenched" theories, i.e. those in which there are internal quark loops, but the 
masses of the valence and loop quarks differ. 

The outline of this paper is as follows. In the following section, we develop the quenched 
chiral Lagrangian for baryons. We then, in Sec. ^|, present the Feynman rules and sketch a 
sample calculation. Section ^ contains our results for baryon mass renormalizations. 

2 Chiral Lagrangian for Quenched Baryons 

We calculate the quenched chiral corrections to the baryon masses using an effective La- 
grangian of "heavy" baryons coupled to the low-lying pseudoscalars. The approach is a 
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synthesis of two techniques which have recently appeared in the literature: 1) Baryons 
are treated as fixed velocity sources, thereby eliminating the baryon mass term from the 
Lagrangian and 2) The theory is extended from QCD to QQCD by constructing a 
Lagrangian symmetric under the graded group SU(3|3)l x SU(3|3)# M. The latter step 
corresponds to the addition of "bosonic quark" degrees of freedom whose internal loops in 
Feynman diagrams cancel those of the ordinary quarks. 



2.1 Review of results for QCD 

We start with the low-energy effective theory for QCD, in which the pseudo-Goldstone bosons 
are coupled to octet and decuplet baryons. The Goldstone fields are grouped in the usual 
3x3 matrix ir(x) (normalized so that 7T%3 = The octet baryons are similarly 

grouped in the 3x3 Dirac field B(x) (with _B 13 normalized to the proton). Finally, the spin- 
3/2 baryons are represented by a Rarita-Schwinger tensor T^^x), symmetric in its three 
flavor indices (with Tin normalized to the A ++ ). Fixed velocity fields are defined by 

B v {x) = ex.Y>{imBfv^x^)B{x) , T£(x) = exp(ims^u Al x /i )T^(x) , (3) 

where v is the velocity, and mg is the mass of the octet baryons in the chiral limit. As shown 
in Ref. @], in the mg — ► oo limit, the Dirac gamma matrix structure of the theory can be 
eliminated in favor of and the spin-operators S%. We will explain the properties of S£ as 
they are needed. 

To construct a chirally invariant Lagrangian, we use fields which have simple chiral trans- 
formation properties. These are the exponentiated pion fields 

E(x) = exp(2z7r(x)//) , £(x) = exp(i7r(x)//) , (4) 

the axial and vector currents 

A, = ~ , V, = i(£c^t + , (5) 

the mass terms 

M ± = e f Me f ± £M£ , M = diag(m M , m d , m s ) , (6) 
and the covariant derivatives of baryon fields 

WB = d^B + [V M , B] , (7) 

2>"2& = + {v)iT«ik + n-") : fv;;., + (v*)*^ . (s) 

Here, as in the following, we have dropped the subscript v on the heavy baryon fields — it is 
always implicitly present. Under SU(3)i x SU(3)# the meson fields transform as 

S -> TSi? f , (9) 
Z^LZrt(x) = U(x)£I$, (10) 

where U(x) is defined implicitly through Eq. ([10f). A^, J^A ± , B v and T) li B v ^ which are octets 
under the diagonal SU(3), all transform like 

B v (x) -> U{x)B v U\x) . (11) 
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Finally, the decuplet field and its covariant derivative both transform as tensors, e.g. 

Tijk ~^ UiiU jm U kn T^ mn . (12) 

Here, and in the following, we do not distinguish between raised and lowered flavor indices. 
It will always be clear which indices transform with U and which with W . 

The most general chirally-symmetric Lagrangian can now be written down as an expan- 
sion in momenta and quark masses. The leading analytic and non-analytic corrections to 
the baryon masses are obtained from the following terms 

C v = if 2 tr(^E9 M Et + 2/i M + ) , (13) 
C B n = itr(BvVB) 

+ 2D tr (BS^A^B}) + 2F ^(BS^A^B]) 

+ 2^b D tr(B{M + , B}) + 2yb F tv(B[M + , B]) + 2/i6 tr(BB) tr(M + ) , (14) 
Cm = -iT{v-V)T v + AMTT V 

+ 2HT U S^A I1 T U + C {T V A V B + BA V T U ) . 

+ cTM + T v - aT v T v ti{M + ) . (15) 

We have followed the notation of Jenkins and Manohar ||, except for the octet mass terms 
in which we follow Bernard et al. ||. The traces involving the octet baryon are over flavor 
indices. The contractions of flavor indices for terms involving decuplet baryons are not 
shown explicitly — we discuss them for the quenched Lagrangian below. Finally, AM is the 
decuplet-octet mass splitting in the chiral limit. 



2.2 Quenched chiral Lagrangian for mesons 



We now consider the quenched theory. Bernard and Golterman have developed a Lagrangian 
framework which provides a consistent means for calculating the physics of the low-lying 
pseudoscalars in QQCD ||. We briefly review this technique, introduce a compatible rep- 
resentation for the baryons, and then construct the quenched baryon Lagrangian analogous 



to that of Eqs. (0>HD above. 

QQCD can be described by the addition of bosonic quark degrees of freedom, <f i5 one for 
each flavor of light quark. They have the same masses, one for one, as the original quarks, 
qt. Integrating over the qi in the functional integral yields a determinant which exactly 
cancels that resulting from the quark integration. This theory is symmetric, at the classical 
level, under the graded group U(3|3)l x U(3|3)#, and this symmetry dictates the form of the 
low-energy effective theory for the pseudo-Goldstone bosons. To construct this theory, one 
replaces the field ir(x) in the definitions (f|) by the field <&(x), given in block matrix form as 



X 7T 



(16) 



Here ix contains the ordinary mesons (qq), vr the mesons composed of bosonic quarks (qq), and 
X and t ne "fermionic mesons" (qq and qq respectively). The mass matrix is extended to 
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M = diag(m u ,md,m s ,m u ,md,m s ), while the definitions of axial and vector currents retain 
the same form (|5|). The transformation properties of the fields are unchanged, except that 
the matrices L, R and U are now elements of Z7(3|3). 

Because of the anomaly, the chiral symmetry of the quantum theory is reduced to the 
semidirect product [SU(3|3)^ x SU(3|3)#] <8> U(l). The reduction in symmetry allows the 
quenched chiral Lagrangian to contain arbitrary functions of the field 

$0 = ^8^=^(^-77'), (17) 

(7/ is the usual SU(3) singlet meson, rj' its ghostly counterpart), since $0 is invariant under 
the quantum, though not the classical, symmetry group. Putting this all together, one arrives 
at the following quenched replacement for Eq. (|T3|) 



AQ) f 



* 4 



2 

W^E^E*) Ki($ ) + 2/i str(.M + ) V 2 (%)] ( I ? 



where str denotes supertrace. The potentials can be chosen to be even functions of $ , an d 
are normalized as Vi($o) — 1 + C(^o)- We w ^ no ^ * n feet? need the higher order terms 
in these potentials. Our notation follows that of Bernard and Golterman ||, except that 
we use m instead of /i (following Ref. ||) and a$ instead of a. We also differ from M in 
choosing the normalization of / such that /„- ~ 93 MeV. 

At this point the development diverges from that in QCD. In QCD, the last two terms in 
lead to wavefunction renormalization and a mass shift for the 77'. The 77' is then heavy, 
and can be integrated out of the theory, giving the usual Lagrangian C n fll3]). In QQCD, 
by contrast, there is a cancelation between diagrams with more than one insertion of either 
a$ or m on an 77' (or 77') propagator || |], pLOfl . Thus the 77' and 77' remain light, and the 
a$ and m\ terms must be included as interactions in C^K This leads to a more singular 
behavior of a number of quantities (e.g. m n and fx) in the chiral limit. Furthermore, the 
new vertices destroy the usual power counting. Higher loop diagrams involving a$ and mo 
are not suppressed by powers of p/A x or m n /A x , where A x w 1 GeV is the chiral cut-off. We 
assume that a$ and m /A x are small, and work only to first order in these parameters. 

Finally, although we have been talking about the 77', it is not, in fact, a mass eigenstate 
unless the quarks are degenerate. Instead, the flavor-neutral eigenstates are those with flavor 
composition uu, old and ss, with squared masses M^ q = 2/im g (q = u, d, s). These replace 
the 7To and 77 of QCD. The flavor non-diagonal mesons are the same in both theories, having 
the form g^-, with Mfj = fi(rrii + rrij). 



2.3 Baryon representations in quenched QCD 

To construct the chiral Lagrangian for baryons in QQCD, we need to generalize the octet and 
decuplet fields B and T M . The corresponding quenched fields will contain additional baryons 
with compositions qqq, qqq and qqq. Even if we restrict the external states to be the usual 
qqq baryons, the extra baryons will appear in loops. Thus we must include these states in our 
effective Lagrangian. Just as when constructing the baryon part of the chiral Lagrangian in 
QCD, all we need to know is how the states transform under the vector subgroup SU(3|3)y 
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(i.e. L = R = U). The transformations under the full group are simply obtained by using 
the position dependent U{x) defined in Eq. fllPp. Thus we need to determine the irreducible 
representations of SU(3|3) which, when restricted to the quark sector, contain only an octet 
or a decuplet of SU(3). 

We begin with the octet, which is the more difficult case since it has mixed symmetry. 
We construct representations using the "quark" field Q = (u, d, s, u, d, s) and its conjugate 
Q. Under SU(3|3)y these transform as fundamental and anti-fundamental representations, 
respectively 

Qi-^U lj Q j and Q.—.Q.U^ (i,j = l,6). (19) 

We now define the tensor spin-1/2 baryon field Bijk{x) to have the same transformation 
properties as the following operator constructed from the generalized quark fields 



TV 



QT a QTQT ~ Qi' a Q]' c Qk b e abc (Ci 5 U • (20) 



Here C = 27270 is the charge conjugation matrix, and a, b and c are color indices. We 
have raised the color and spinor indices on the fields for the sake of clarity, and suppressed 
the common position argument of all fields. The transformation of Bij k under SU(3\3)v is 
defined through the r.h.s of Eq. (|20|) — the Q's are first rotated "inside" the operator, and 
then the U's are moved to the outside, giving rise to a grading factor because the off-diagonal 
3x3 blocks of U are Grassman variables. The result is 

B] jk — U^WW^U^UwUmIQw . (21) 

Here we are using the following notation in the grading factor: if the index on the field is 
"anticommuting" (in the range 1 — 3) then the index equals 1 in the grading factor; if the 
field index is "commuting" (in the range 4 — 6) the corresponding index equals in the 
grading factor. 

This construction of B automatically yields a representation of SU(3\3)y, since it is 
written in terms of Q's. The Dirac and color structure of the operator enforce the constraints 
coming from the fact that B creates spin-1/2 baryons which are color singlets. The second 
term on the r.h.s. of Eq. ( p0|) is required in order to make the representation irreducible. To 
see this, consider the operator which result when the indices % — k are restricted to lie in the 
range 1 — 3: 

B ijk \ R = B! jk + B! kj (22) 

(R indicating restriction), where the usual quark baryon operator is defined to transform as 
in 

&ijk ~ [Qi^Qj' b Qk' C \ £abc(C~/ 5 ) a p . (23) 

B q is anti-symmetric under % <-> j, which allows the first two quark indices to be combined 
into an antiquark index 

r., 

ijk 

Now the problem is clear — the operator B k creates both SU (3) octets and singlets, but we 
only want the former. The singlet is obtained by contracting B q with Eij k . It can be cancelled 
by symmetrizing the last two indices of B q , as is done in Eq. (^). The octet part is, up to 
an overall constant, just the standard 3x3 matrix used to represent baryons in the QCD 



B% k oc eij „Bi . (24) 
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chiral Lagrangian. Using the normalizations discussed below, the explicit relation between 
B q and this field is 

^ijk R = 7e \ £ ijk' B k + £ ikk'Bj ) . (25) 

We have convinced ourselves that B transforms in an irreducible representation of SU (3|3)y. 
It has dimension 70, and decomposes under SU(3) as follows: an 8 each of qqq and qqq states, 
and al + 8 + 8 + 10 each of qqq and qqq states. It satisfies the symmetry properties (easily 
determined from the above definitions) 

Bijk ( ) Bikj , 

= B ijk + (-T j+1 B jik + (-) ij+jk+ki+1 B kji . (26) 

These relations show why B cannot be reduced to a two index form in QQCD — it has no 
grading independent symmetry properties under the interchange of two indices. 

The need for a three index tensor makes sense also from another point of view. The major 
aim of the formalism we are constructing is to allow identification of contributions to mass 
renormalization which contain internal quark loops. To do this we need to be able to follow 
the flow of flavor through the Feynman diagrams. This is not straightforward using the 
standard 2-index form for B, since one of the indices carries the flavor of two quarks. Using 
the 3- index field B, by contrast, there is a one-to-one correspondence between the terms 
which comprise a given diagram in QChPT and the flow of quark flavor. In practice what 
happens is that all diagrams in which the quark flow contains internal loops are cancelled 
by the corresponding diagrams with q in the loops. 

The representation containing the spin-3/2 baryons is simpler to construct. The decuplet 
baryons are already represented by a symmetric three-index tensor in QCD. One only needs 
to extend the range of the indices, and apply the appropriate symmetrization. We find that 



a,ijk 



QTQfQY + Qi' b QTQk a + QTQTQk h ] e^JC-f)^ . (27) 



This transforms under SU(3|3)y exactly as does Bijk, Eq. (]2lj) . It has the following symmetry 
properties (dropping Lorentz and Dirac indices for clarity) 

r m = (-) ij+1 T jik = (-y k+1 r tkJ . (28) 

These imply that the representation is 38 dimensional, containing a 10 of qqq states, a 10 
and an 8 of qqcfs, an 8 and a 1 of qqtfs., and a lone 1 qqq. When the indices are restricted 
to the range 1 — 3, the QQCD tensor is proportional to the decuplet tensor used in QCD, 
and we choose normalizations so that 

Tijk n T-ijk ■ (29) 
H 



2.4 Quenched chiral Lagrangian for baryons 

To construct the quenched generalization of Eqs. ([14]) and ( |15D we need quantities which 
are invariant under SU(3|3)^ x SU(3|3)#. These we construct from the fields B and T, their 
covariant derivatives, and their conjugates. The covariant derivatives of both fields take the 
same form, exemplified by 

V»B ijk = d»B ijk + (vn«S m + {-) w) {V% r B ifk + (-f+^ k + k ')(y») kk ,B l3k> , (30) 
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where V* is the vector-current defined by the quenched generalization of Eq. (|5]). The 
grading factors arise because, in order to compensate for the position dependence of the U's, 
the V's must act "inside" B — i.e. as if they were coupled directly to the generalized quark 
fields in Eq. |2D|. These covariant derivatives transform under SU(3|3)l x SU(3|3)# in the 
same way as the fields B and T. The conjugate baryon fields are defined by 



B 



Vfc Vj Vi Vfc V 7 Vj 



a/3 



T 



kji, a 



Ql Q7 Qi + Qk Qj Qi + Qk Qi Ql 



and both transform in the same way, e.g. 



B 



\ l ) °i'j'k' u k'k u j'j u i'i ■ 



(31) 
(32) 

(33) 



Using the group transformation properties, it is straightforward to catalogue all the 
bilinear invariants. We use the compact notation 



ijk,/3 j 



(STB) = K^a P B m 
(BTAB) 

(BTBA) = % x T a p A kW B ijV , p x (- 
(7**I\4%) 



B a kji rjA u ,B 
>kji 

Trf 1 -pa rj-fi 
kji,a fi,ijk ' 

1 kji,a L 0-™-W 1 ^L,i'jk > 

~j2 a -pa A[l rrfi 
'-'kji 1 P-n-ii' I fj,,i'jk ■ 



Ji+j)(k+k>) 



(34) 
(35) 
(36) 
(37) 
(38) 
(39) 



There is also the conjugate of the last quantity, i.e. (T^YA^B). In the above constructions, 
A is an operator transforming like the axial current A^, and T is an arbitrary Dirac matrix. 
In fact, the change to fixed-velocity fields simplifies the Dirac structure as for QCD, and in 
practice only the spin operator enters in place of T. Various possible terms are absent 
from the list above because they can be rewritten in the above form using the symmetries 
of B and T. For example, a (BAB) bilinear with A coupling to the second index of Bijk 
is redundant due to the symmetry B^k = (—\y k+1 Bikj- Other terms simply vanish, e.g. 
constructions involving BkjiTtjk- 

We can now write down the relevant part of the quenched Lagrangian for baryons. It 
consists of 



C 



= i(BvVB) 

+ 2a (BS^BA^) + 2/5 (BS^A^B) 
+ a M (BBM + ) + (3 M (BM + B) - 



- 2 7 (BS^B)str(A IJ ) 
a (BB)str(M + ) 



(40) 



and 



-%(j\v-V)%) + AM(T U %) 



+ 2H(T U S' 1 A U T P 



2 U 



(T V A U B) + (BA V T V )\ + 2i{T v S^T v )str{A v ) (41) 



+ c(T u M + %) - a (T u %) str(M+) 
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Each term can be multiplied by an arbitrary even function of $o ; but the higher order 
vertices these functions produce do not contribute at the order we are working. There are 
also terms proportional to (£>£>)$q an d (BB)d^&od^&o, but these do not contribute to mass 
renormalization at one-loop. In quark-flow language, these are double hairpin vertices and 
lead to closed quark loops 0. Finally we note that str(A/) is proportional to d u & , so 
possible additional terms involving the latter are not independent. 

The quenched Lagrangian ([4T)l)-(|4~l|) looks very similar in form to that for QCD, Eqs. 
(pT4|)- (p~5|) . To make the correspondence precise we need to pick out the parts of and 
£r$ containing only qqq fields. This is straightforward for the terms involving the spin-3/2 
field, since both quenched and full fields have three indices, so the structure of the allowed 
contractions is the same. For example, using Eq. (p9[) , one finds that 

(T U S»A,T V ) 



T V S»A^T U . (42) 



Thus the coefficients AM, c and 7i play the same role in the quenched Lagrangian as they 
do in the full theory. It is important to realize, however, that there is no reason for the 
coefficients to have the same values in the two theories. Despite this, we use the same 
symbols so that the physical significance of each term can be more easily recognized. 

Terms involving the spin-1/2 field require more work to interpret. We need to convert 
from the three index tensor B to the matrix B. Using Eq. fl2"5|) we find 



(BB) r = tr(BB) , (43) 

(BBA)\ r = ltx(BAB) + ltr(BB)tx(A) - \ti(BBA) , (44) 

(BAB) | = -\tr(BAB) + |tr(S J B)tr(Aj - ftr^BAj , (45) 

(EB)str(Aj = tr(5S)tr(A) . (46) 



(T V A U B) r = -«1\T V A V B = -^lT v ijk A u>iil B jf s ilfk . (47) 

The first of these relations is actually the condition which sets the normalization in (P5]). 
From Eq. (f|7|) we see that the coefficient C has the same significance in quenched and full 
QCD. As for the other coefficients, the quenched Lagrangian is equal to that for QCD, Eq. 
(|14|), if we make the identifications 

a = 2(iD + F), = (-§£> + F), 1 = 2{D-F), 

i ( 48 ) 
a M = 4(^b D + b F )fj,, (3 M = 2(-|&d + b F )H, a = 2(b + b D - b F )fi . 

With the exception of the result for 7, we use these relations to reexpress our quenched 
results in terms of D, F, bp, bp and 60, rather than a, (3, cum, Pm and a. This allows a 
more direct comparison with results from chiral perturbation theory in QCD. We reiterate, 
however, that the values of the coefficients will be different in the two theories. 

With all these correspondences in hand, we can now see the most important difference 
between the Lagrangians in QCD and QQCD, namely that the latter has two additional 
coefficients, 7 and 7'. Eq. ( |48D shows that 7 is non-zero in QCD, but that it is not independent 
of D and F. In QQCD, by contrast, 7 is an independent parameter. Both the new terms 
involve a baryon bilinear coupled to sti{A u ) oc d u Q , and thus represent couplings of the rj' 



9 



and rf to baryons. These are independent couplings because there is no symmetry connecting 
the couplings of SU(3) octet and singlet mesons. In the quark-flow language, these are 
"hairpin" vertices — the q and q (or q and q) in the rj' (or fj') annihilate. Such couplings are 
present in QCD, but are absent from the chiral Lagrangian because the rf is heavy and can 
be integrated out. 



3 Flavor Feynman rules and a sample calculation 

Calculations of Feynman graphs separate into "Dirac" and flavor parts. The former are 
standard, using the Feynman rules and methodology worked out in Refs. |J. We do not 
give further details here. The flavor part of the calculation is new, and we have developed a 
diagrammatic method for carrying it out. This method not only simplifies the calculations, 
but also allows one to trace the underlying quark flows. In this section we explain our 
method, and sketch a sample calculation. 



3.1 Feynman rules 

The quenched baryon Lagrangian is written in terms of the three index tensor fields and 
%jki each of which has 216 components. The symmetry relations [Eqs. ( pB"D and (|28|)1 reduce 
the number of independent fields to 70 and 38, respectively. One could proceed by expressing 
the Lagrangian in terms of these independent fields. By construction the baryon propagators 
would then be diagonal and normalized. The vertices, on the other hand, would involve 
complicated flavor factors. In practice, we find it more useful to reexpress the Feynman 
rules in terms of all of the components of B and T, and impose the relations between the 
components only on external lines. As we show explicitly below, the symmetry relations are 
automatically maintained on internal lines by the structure of the baryon propagators. 

We are interested in diagrams in which the external baryons are qqq states. Thus, we 
need to display the independent fields explicitly when the indices of B and T lie in the range 
1—3. Using eqs. (E5T), the octet baryons are 



p = 


= bu , 


n = 


= 621 , 


£+ 


= b%3 , 




= h+) 


-0 


= 631 , 




= b 32 


A 


= h-) ■ 







'23 , 



(49) 



where we have used the symmetry relations to define 

b tj = VgBuj = VeBiji = -yi (* ^ j), (so) 

6 (+) = V3(B 123 + B 23 i) = -^3B 3 i2 = V3(B 132 + B 213 ) = -V3B 321 , (51) 

&(-) = (B123 - £231) = (B132 - B 213 ) . (52) 

For the decuplet baryons the relation to physical fields is just as for the usual tensor field, 
for example 

A++ = Tin , E*> + = T 113 = T 13 i = T 3U , 
Q~ = %33 , £*' = T 133 = a/3 T 3V3 = a/3 T 331 . 
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We now present the Feynman rules, beginning with the "pion" propagator 



$ n 

m 



- - # = 



P 



(54) 



Here, lines joining indices denote Kronecker deltas, so that the flavor factor is S ni 5 m j. The 
form of the grading factor follows from the the properties of the supertrace. 

To derive the baryon propagators, we write each component of B or T in terms of the 
independent fields, and then use the fact that the propagator is diagonal for these fields. A 
useful trick for obtaining the relative factors, and particularly the grading factors, is to think 
in terms of the underlying Q fields. For example, for the spin-1/2 propagator, we rewrite 
the field as 

Bijk ~ [QiQjQk - (-l) jk QiQkQj} x spin x color, (55) 

insert this relation and its conjugate into the propagator {Bi mn ~B~ l '' k ), and then compute the 
ordinary Wick contractions of the Q vectors. Reordering the Q's to form singlet combinations 
results in a common spin and color factor multiplied by a constant which includes the grading 
factors. In this way we obtain 



B, 



Imn 



&kji 



(1/2) 



lmn,ijk 



Imn ^ 



kji 



T (3/2) 



1 77(1/2) 

{vk)+ie lmn > ijk ' 



l m n 



I m n 



I m n 



2(_)i* 



(56) 
(57) 



i j k i j k i j k 

Imn Imn Imn 

i j k i j k i j k 



ft 



(3/2) 



lmn,ijk 



v-k) + ie 

Imn Imn Imn 

- n jk | X - n ij 

i j k i j k i j k 



(58) 
(59) 



+ (") 



ij+ik 



Imn Imn 

_|_ ^_yk+jk _ ^_yj+ik+jk 

i j k i j k 



I m n 
i j k 



Here P^ u = v^v u — g^ u — (A/3)S fM S 1 ' is the spin-3/2 projection matrix, and the are flavor 
projectors, with j labeling the spin. We discuss the properties of the fts below. 
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We obtain the vertices in a similar fashion. The baryon-pion vertices are 



B 



nml 



Y 



k-S 



2a(- 



Ai+j)(p+q) 



I m n p q 



i j k 



q p I m n 



I m n p q 



+ 2(3 



+ 2 7 



i j k 



i j k 



(60) 




(61) 



Y 



9fj.u k ■ S 



2H 



2i 



(62) 



where k is the incoming pion momentum, and fi, v are the Lorentz indices of the Rarita- 
Schwinger fields. We have shown all flavor indices explicitly for the first vertex in order to 
illustrate the conventions that we use in subsequent vertices. The two-point vertices are 
more straightforward 



X- 



i(ak 



rrin 



(63) 



2a r 



+ 2{3 n 



(64) 



lyu = ig^AM 



(65) 



ig^ y 2c 



(66) 



where we have adopted the notation that a dot on a line with flavor j indicates a factor of 
rrij. For example, the spin 1/2 mass term is 

% S u 8 jm S kn {2a m m k + 2j3 m mi) . 

Note that potential mass terms involving a and Zf vanish because str(M) = 0. In the 
"hairpin" vertex, (|6~3"|), we have included a factor of 2 resulting from the possible contractions. 
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(a) (b) 



Figure 1: Examples of diagrams contributing to spin 1/2 baryon mass and wavefunction 
renormalization. 

We also need vertices in which two pions emanate from a baryon mass term, but for the 
sake of brevity we do not give the results here. They are straightforward extensions of the 
vertex in Eq. (|66|) . 

3.2 A sample calculation 

To illustrate our method, we sketch the computation of mass renormalization of the spin-1/2 
baryons resulting from the diagrams of Fig. ([!]). We mainly focus on diagram (a), and in 
particular the part proportional to (3 2 . The Feynman integral is 

r d 4 k i 2 (k^)(-k u ) 
v J (27r)*(vk)(k 2 -M 2 ) ' 1 ' 

and multiplies uS^S^u, u being the spinor of the external state. For the moment we denote 
the meson mass in the loop generically as M T . The finite part of the integral is 1^ = 
—ig Mi/ M 7r /(247r), allowing us to eliminate the spin vectors using uS 2 u = —3/4. The mass 
renormalization, Am, is thus given by 

M 3 

- i Am = t(3 2 —^ J (68) 

aside from flavor factors. 

The pion in the loop can, for non-degenerate quarks, have one of six masses: M 2 ^ = fi(rrii+ 
rrij) where i < j < 3. This is true even if the pion contains bosonic-(anti) quarks, since they 
are degenerate with the corresponding quarks. In our calculation of spin-3/2 baryon mass 
shifts we have worked with completely non-degenerate quarks. For the spin-1/2 baryons, 
on the other hand, we have considered only the isospin symmetric limit, m = m u = m^. 
This reduces to three the values of the meson masses: M 2 = 2/im, M\ = /i(m + m s ) and 
M S s = 2fim s . The full expression for mass renormalization is then 

M 3 

I=1T,K,SS ■> 

where B specifies the external spin-1/2 baryon, and cb,i is the flavor factor. 

The flavor factors are obtained by projecting the contribution from the flavor part of the 
Feynman rules onto the particular flavor of the external baryon. Explicitly, we can write 

c B j = Y.^DM, (70) 

a.b 
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where the indices a and b run over all 216 values of the three flavor indices in the dependent 
field basis, D 1 ^ is the flavor part of the Feynman graph in this basis, and B is one of the 
independent baryons, whose "wavefuction" is i\) B . The wavefunctions can be read off from 
Eqs. (^Uj)-(^). The flavor "charged" octet baryons are exemplified by 

^112 V>121 ^211 ^332 ^323 ^233 

V=\f\(^ 1> -2), r=^\{\ 1, -2), (71) 
where we show only the non-zero elements. The flavor-neutral octets have wavefuctions 

V>123 ^231 ^312 ^132 ^213 ^321 

^ = \fi{ h 1, -2, 1, 1, -2), (72) 

^123 ^231 ^132 ^213 

^ A = \ ( 1, -1, 1, -l). (73) 

In other calculations, we also need the wavefuctions for the decuplet baryons, ipl- These are 
exemplified by 

1p lu _ ^113 ^131 ^311 

^ A++ = ( 1 ) , ^ + = V?( L L 1 ) > ( 74 ) 

^123 ^231 ^312 ^132 ^213 ^321 
= V / J( 1, 1, 1, 1, 1, 1 ) . (75) 

Explicit forms for the wavefunctions involving bosonic-quarks are not needed, since these do 
not appear as external states in quenched lattice calculations. 

The wavefuctions satisfy various useful properties. Those for the spin-1/2 and spin-3/2 
baryons are separately othonormal 

E «' = { J I ^ % . E Hi? = { \ I J T r ■ (76) 

This ensures that the kinetic terms are correctly normalized when written in terms of in- 
dependent fields. The wavefuctions are eigenvectors of the flavor propagators with unit 
eigenvalues 

All other eigenvectors have zero eigenvalues, corresponding to the fact that the flavor prop- 
agators are projection matrices, i.e. [JF^] 2 = T^' for j = 1/2 and 3/2. Thus the flavor 
propagators can be written as 

^ 1/a) = EtfcW, ^ /2) = £€Vr (78) 

B T 

The sums here are over the independent baryons of the appropriate spin (70 and 38 respec- 
tively), including baryons containing bosonic-quarks. These results show that the baryon 
propagators in Feynman diagrams project onto the independent states, as claimed above. 
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Returning to the sample calculation, we next compute the "matrix elements" D ! ab using 
the flavor Feynman rules. We begin with the elements needed if the external baryon is a 
proton, i.e. a,b— 112, 121,211. A representative element is 



12,112 



E 



1={1, 2,4,5} 



112 



+ 



+ 



1 1 2 



112 



(79) 



1 

6 



E (-i) m 

£={1,1,0,0} 



1 

12 



1 

12 



1 

6' 



D K 

-^112,112 



E 

'={3,6} 
'i={l,0> 



1 1 2 
.1)1+1 




Dtkiu = (81) 

In these equations we are using the diagrammatic notation in a slightly different way from 
above — the diagrams denote both the way in which the flavor indices are contracted and the 
corresponding factors from the propagators and vertices. The remaining elements can be 
evaluated similarly, with the results 



£>* = -- I 1 1 I , D K = 0, D ss = 0, (82) 
6 \0 2) 

for the (112, 121,211) block. Sandwiching this between the proton wavefunctions gives the 
flavor factors c p j = —1/3, 0, 0, for I = n,K, ss. These can be inserted in Eq. (]69|) to obtain 
the proton mass renormalization. 

It is worthwhile contrasting the calculation thus far with that in QCD. To make the 
comparison, it is better to use a three index field to represent spin-1/2 baryons in QCD, 
for then one can develop the calculation diagrammatically as above. The major difference 
from QQCD is that the indices run from 1-3 instead of 1-6. (One must also project against 
the 7]' in the meson propagator in QCD, but this is a small effect numerically.) Thus the 
main difference between calculations in the two theories is that quark-flow diagrams involving 
internal loops cancel in QQCD, but do not in QCD. This is, qualitatively, exactly as expected. 
What we are able to do here is make a quantitative calculation of the effect of this cancelation. 

A striking feature of the calculation for baryons is that there are quark-flow diagrams 
which contribute in QCD which survive the cancelation between q and q loops, e.g. the last 
two diagrams contributing to D\ V1XV1 . This is in contrast to mesonic quantities (e.g. 
and /„-), where the cancelation of quark-loops removes all the QCD contributions ||, [|. The 
surviving diagrams for baryons must, however, contain mesons composed of quarks having 



15 



the flavors of the valence quarks of the external baryon. Thus, since the proton does not 
contain a valence strange quark, there can be no contributions from kaons or ss mesons in 
the loop. If the external particle contains a strange quark, however, kaon contributions are 
present, e.g. 



r)7r 

^113,113 



-^131.113 



1 1 3 




1 

12' 
1 

12 



-^113,113 



-^131,113 




113 




1 

12' 
1 

12 ' 



1 1 3 



113 



Here we have not shown the quark-loop cancellations explicitly. 

Proceeding as above, we calculate all the relevant parts of D 1 . Contracting these with 
the wavef unctions, we obtain the flavor factors 



18 c B , 





7T 


K 


ss 


p 







°\ 


s 


1 


5 










5 


1 


A 


\3 


3 


o J 



J3) 



When inserted in Eq. (^9|) this gives the result for the contribution proportional to (3 coming 
from Fig. (TTja) . 

The calculation of the other contributions proceeds similarly. The diagrams for the a 2 
terms are the same as those just discussed, except that the loops come off the right hand 
side of the baryon instead of the left. This changes the flavor factors. The diagrams for the 
a(3 contributions are different, and are shown in Fig. [2]. Again, some of the diagrams involve 
internal quark loops and cancel in the quenched approximation, whereas others remain. 
Diagrams for the ja contributions are shown in Fig. |3|. These are special to the quenched 
approximation because only flavor singlet pions can appear in the loop. The same is true 
of contributions proportional to 7/?. Finally, terms proportional to 7 2 vanish because they 
involve internal quark loops. 

Additional types of quark flow occur in the evaluation of Fig. (|l]b). This diagram is 



not present in QCD because only flavor singlet pions couple to the 



(or a$) vertex. 



Examples of the corresponding quark-flow diagrams, for contributions proportional to (3 Wg, 
are shown in Fig. There are similar diagrams for the a 2 and a(3 contributions. It is these 
type of diagrams which are the sole quenched contribution to mesonic quantities. Finally, 
contributions proportional to 7 involve quark loops and cancel. 



3.3 Quenched oddities 

We close this section by noting some peculiar features of the workings of chiral perturba- 
tion theory for QQCD. We first show why care must be taken if one attempts to calculate 
quenched results without the benefit of the graded formalism. This might be done, for ex- 
ample, by making a quark-flow correspondence using the QCD chiral Lagrangian. Such an 
approach was used successfully for mesons in Ref. |§. 
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Figure 2: Quark-flow diagrams contributing to the a/3 terms in mass renormalization. Arrows 
have been added to help clarify quark flows. The boldface line can be both a q and a q; 
other lines represent only quarks. 




Figure 3: Quark-flow diagrams contributing to the 7a terms in mass renormalization. 




Figure 4: Quark-flow diagrams contributing to the (3 2 ml terms in mass renormalization. 
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(a) 



(b) 



Figure 5: Diagrams contributing to the A self-energy which can lead to A decay. 



Consider D 211 2 n , which can be split up as follows 



-^211, 



211 



+ 



2 11 



2 11 




2 11 




+ 



J4) 



2 11 2 11 

The correct approach is to drop all diagrams containing quark loops, since these cancel in 
pairs. This leaves the last four diagrams on the first line, which give a non-zero result, —1/3. 

One might have been tempted, however, to argue as follows. In each of these four 
diagrams, the baryon propagating in the loop is 8m. There is, however, no spin-1/2 baryon 
consisting of three up quarks — indeed, the symmetry relations set £>m = 0. Thus -Dfn 211 
should vanish in the quenched approximation. 

This argument is wrong because, as we have seen, the baryon propagator automatically 
projects against B\\\. What happens is that the six diagrams on the first line (the only ones 
containing £>m) cancel: their sum is proportional to (2 + 2 — 1 — 1 — 1 — 1) =0. In the second 
line, the contributions from I = 2 and I = 5 cancel as usual, leaving the entire contribution 
to come from the diagram with I = 4. This consists entirely of diagrams containing an 
internal u loop! Thus, if one insists on removing £>m, then to obtain the correct result in 
the quenched theory one must include bosonic quarks. Clearly the simplest approach is to 
let the theory take care of doing the projections itself. 

A related peculiarity concerns the decay of decuplet baryons in QQCD. In QCD they 
decay through the strong interactions, e.g. A ++ — > n + p, due to the vertex proportional to C. 
One could study this numerically by calculating the A two-point function in Euclidean space. 
At large Euclidean times it would be dominated by the n + p intermediate state, resulting 
from the cut in Fig. (0a) .0 The issue is whether the same is true in QQCD. The following 
argument suggests that it is not. The A ++ has quark composition uuu, and, since there are 
no quark loops in QQCD, it can decay only into a uu meson and a uuu spin-1/2 baryon. 
But there is no uuu octet baryon, and so the decay cannot occur. By isospin symmetry, all 
of the A's must be stable in QQCD. 

This argument is wrong for the same reason that -DJn 211 does not vanish. One can only 
impose the requirement that the Bm does not contribute if one works with the bosonic-quark 

2 Strictly speaking, this is only true in large enough volumes. The decay amplitude is p-wave, and so 
vanishes at threshold, requiring the intermediate pion to have non-zero momentum. This increases the 
threshold for the decay on small lattices. 
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formalism in which all symmetries are manifest. But then there are additional contributions 
to the self-energy which do lead the A to decay. If one instead simply throws away all quark 
loops, then the uuu baryon does contribute in intermediate states, and again the A decays. 
Another way of understanding this is to note that quark loops are necessary in order to 
implement the Pauli exclusion principle on internal propagators. 

Thus we claim that one can study the decay of the A even in QQCD. This is such a 
counter-intuitive result that it is worthwhile understanding in more detail. Consider the 
quark flow diagrams contributing to the C 2 term in Fig. (|]a) 



"^111,1 



11 



+ 



+ 




+ 



1 1 1 



111 



111 



111 



111 



111 



+ E 



Z={2,4,5} 



+ (J 



55) 



hi 



ill 



The quark-flows are the same as those contributing to the (3 2 part of D% n 211 — the only 
difference is the external flavors. The story is now the same as above. The diagrams with 
quark loops cancel, leaving a non-zero contribution from the last four diagrams on the first 
line — all of which involve Bm. Alternatively, the Pauli exclusion principle requires that 
the diagrams on the first line cancel, leaving a non-zero contribution from the diagrams of 
the second line with £ = 4. In this view, the quenched decay is actually A ++ — > (uu) + 
(uuu). In QCD, by contrast, the diagrams on the first line still cancel, but the non- vanishing 
contribution is from i = 2 on the second line. Thus the QCD and QQCD decay amplitudes 
have the same magnitude but opposite sign (assuming m u = m d ). 

Actually this is not the whole story in QQCD. There are contributions from Fig. (|5]a) 
proportional to £7', and there are contributions from Fig. (|5|b) proportional to C 2 m\ and 
C 2 a%. Both give rise to a pion-nucleon cut in the A propagator. Thus the A will have a 
different width in QQCD than in QCD. Nevertheless, it will be interesting to use QQCD 
as a testing ground for methods to study an unstable A. This requires, of course, that 
wa — m p > m n , a condition which is not close to being satisfied in most present simulations. 
One must use lighter quarks (to increase ?tia — rn p and to decrease mj and larger lattices 
(to decrease the cost of having non-zero momentum pions) for the decay to be kinematically 
allowed. 



4 Results 

Using the methods explained above, we have calculated mass renormalization of octet and 
decuplet baryons in QQCD resulting from the diagrams shown in Fig. |6|. These diagrams 
lead to the dominant corrections in a limit which we explain below. The calculation of the 
flavor factors is quite laborious. In order to check our results, we have done two independent 
calculations, one by hand and the other using Mathematica. We have also tested our method 
by repeating the calculation for QCD, and comparing to existing results. 

We present our results separately for each subset of graphs, i.e. (a) — (e) in Fig. We 
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(a) 



(b) 



(c) 



/ \ / \ 



^x^ ^x^ 



(d) 





(e) 




-x. 



x* 



Figure 6: Diagrams included in this calculation. Those in which a loop is attached to an 
external line [e.g. the first and third diagrams in (d)] represent mass shifts coming from 
wavefunction renormalization. 



use the notation 



M hary (m u , m d , m s ) = M + M b%( m u, m d , m s ) 

G=a,b,c,d,e 



36) 



Here M is the octet baryon mass in the chiral limit. For the octet baryons we have done 
the calculation assuming m u = rrid, while for the decuplets we have allowed all three quarks 
to have different masses. 

A considerable simplification results from the fact that, in QQCD, the mass of a given 
baryon depends only on the masses of the quarks that it contains. This means that we need 
only give results for the £ and A octets, and the (which has composition uds) decuplet. 
Mass shifts for the other baryons can be found simply by changing the quark masses, e.g. 



M s = M E (m n ,m s 



M p = Mx(m u ,m u ) , 



M a = M s (m s ,mJ 



M s *,o = M^*(m u ,m dl m s ) , M£ + = M s *(m u , m u , m u ) , M n = M s *(m s , m s , m s ) . 



An important check on our results is that all octet baryons are degenerate if the quarks are 
degenerate, which implies Ms(m u ,m u ) = M^{m Ul m u ). 

The diagrams of Fig, (||a) give contributions proportional to M^ 



M s a) = -Ab F M 2 uu + 2{b D - b F )M t 



2 

ss i 



(88) 
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M 
M 



(a) 



4(26 D /3 - bp)M 2 u - 2{b D /3 + 6 F )M S 2 
AM + c'(Mi + M 2 dd + M S 2 J , 



(89) 
(90) 



where d = c/(3fi). These results require some explanation. 

1. The diagrams with filled triangles give contributions proportional to quark masses. 
These corrections are the same as in QCD, except that "sigma-term" contributions 
proportional to a or a are absent in QQCD. To make the comparison with QCD more 
direct, we have rewritten the octet results using Eq. (pE8|) to express ozm and (3m in 
terms of bp> and bp. 

2. The decuplet masses are also shifted by AM, the octet-decuplet mass splitting in the 
chiral limit. 

3. The diagrams involving meson loops with the hairpin vertex give corrections propor- 
tional to m 9 mg log(m 2 ), which are quenched artifacts. For small enough m n , these can 
become larger than the analytic corrections proportional to m q . They are of the same 
form as the leading chiral correction to the pion mass H, ffl 



M 



2fim q 



M 2 
x J 



(91) 



where^l 5 = m 2 ,/ '48n 2 f 2 . It turns out that the corrections to baryon masses from these 
loop graphs are of exactly the form to allow them to be absorbed by replacing quark 
masses with the corresponding 1-loop corrected meson squared masses of Eq. (pT|). 
This is convenient, since it is more straightforward to extract meson masses than 
quark masses from simulations. 

The diagrams of Fig. (^p) are present in both QCD and QQCD. As shown in the example 
of Sec. (|^), these give rise to non-analytic terms proportional to M^. We find 



M 



(6) 



M, 



(b) 



M 



(«0 



where 



:p 5 



2F 2 - 4F 7 



±C 2 W 



+ (ID 2 - ADF + 2F 2 - IC 2 )V US + 2(D - F)jV sl 



id 



DF + 2F 2 + 4(2D/3 - F) 7 



±C 2 

3 U 



+ 



f D 2 - IDF 



2F 2 - \C 2 



V us -[2(D/3 + F) 7 ] V s; 



[lC 2 - l £H 2 )[V ud + V ds + V m 

- §ni[v uu + v dd + v ss ] , 



v u . 



M 3 

us 



etc. 



16/% 

Note that we have rewritten a and (3 in terms of D and F using Eq. 



(92) 

(93) 
(94) 

(95) 



3 Strictly speaking the meson mass in the logarithm should be the lowest order approximation (2/im g ), 
but it can be replaced by the actual mass to the order we are working. 
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These results for show that the M% terms are present in QQCD. Some are propor- 
tional to D and F, i.e. of the same general form as those in QCD, while others, proportional 
to 7 and 7' are quenched artifacts. To give some idea of the effect of quenching, we quote 
the QCD result for A § 



M| fe) (QCD)= 2D 2 + C 2 V v -l {D 2 + 9F 2 )+C 



K 



(96) 



There is no obvious correlation between quenched and full results. 

The diagrams of Fig. (^c) are obtained from those just considered by inserting an m 2 , or 
a$ vertex on the meson line. We find 



M. 



M, 



M. 



AF 2 W(u, u) - A(D - F)F W(u, s) + (D - F) 2 W(s, s) 



+ \C 2 [W(u, u) - 2W(u, s) + W(s, s)] 



(4D/3-2F) 2 W(u, u) - 2(4D/3-2F)( J D/3 + F) W(u, s) 
+ (D/3+F) 2 W(s,s) 



(97) 
(98) 



^H 2 [W(u, u) + W(d, d) + W(s, s) + 2W(u, d) + 2W(d, s) + 2W(s, u)} 

+ ±C 2 [W(u, u) + W(d, d) + W(s, s) - W{u, d) - W(d, s) - W(s } u)] , (99) 



where 



and 



W{a,b) 



nr. 



Ml) 



16/ 2 tt(M ( 



aa 



Ml) 



W(a, a) 



limW(a,b) 



m 2 M a 



(100) 



(101) 



16/ 2 7T 

These contributions are pure quenched artifacts. The a$ term is of roughly the same form as 
the "conventional" loop corrections proportional to M%. The m 2 , term, however, is enhanced 
by rn^/M 2 in the chiral limit. It turns out, however, that the influence of these terms on the 
curves of Mb ary versus M 2 q is not significant until quark masses substantially smaller than 
those used in present simulations. We discuss this in detail in the companion paper. 

It is simple to understand the origin of the enhanced loop contribution. The second 
propagator in the loop makes the diagrams more infrared divergent than those in Fig. (Qb), 
and thus more sensitive to the IR cut-off. The contribution proportional to m 2 , can be 
obtained, up to overall factors, by acting with ml(d/dM 2 q ) on the results from Fig. (|§b). 
This replaces one power of M 2 q with m 2 ,. This is another example of the peculiar behavior 
of quenched quantities in the chiral limit, previously seen in M 2 [see Eq. (fJT|)], f K , f B and 
various matrix elements [§, ^, 0, [TT|] . 

The diagrams of Figs. (f|i) and (e) yield further quenched artifacts which are proportional 
to m q m 2 ] \og(M 2 q ). Without the hairpin vertex, these diagrams give corrections proportional 
to m q M 2 log (M 2 ), which are suppressed in the chiral limit. The insertion of the hairpin 
vertex makes the diagrams more IR singular, replacing M 2 with tbq. In quoting our results 
we have converted factors of m q into M qq using the lowest order formula. This is consistent 
since, as we explain below, we are dropping terms of 0(t?7-q). 

We find that 



Ms 



(d) 



M 



id) 



(d) 



0. 



(102) 
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The cancelation between the two diagrams is related to the simple flavor structure of graphs 
involving the m\ vertex — as illustrated in Fig. (f|). There is no such cancelation between the 
diagrams in (e), and we obtain 



-4b F Mi + 2{b D - b F ) M 2 SS - c'(2 M 2 UU + M 2 



x 



[X(u,u)-2X(u,s) + X(s,s)], (103) 
Mi e) = 0, (104) 
M§ = ±C 2 [d + 2(b f - b D /3)} (Mi + M 2 dd + Ml 



x 

[X(u, u) + X(d, d) + X(s, s) - X(u, d) - X{d, s) - X(s, u)] 
- ^C 2 MMi - M 2 dd ) [X(u, u) - X(d, d) - 2 X(u, s) + 2X(d, s)] 
-±C 2 b D (M 2 uu + M 2 dd -2Ml)x 
[X(u, u) + X(d, d) -2X(s,s)- 4I(m, d) + 2 X(u, s) + 2 X(d, s)} . (105) 



Here 



and 



A ( a -") = • (106) 

X M =^X(a, b) = <^ + 1 l (107) 

The last two terms in M-£i are symmetric under permutations, despite appearances. 

We can now explain why we consider only the diagrams of Fig. @. We are essentially 
carrying out a standard chiral expansion in M 2 q , supplemented by an expansion in ml/3. The 
next chiral corrections are ~ M qq log M qq , and thus suppressed by one power of M qq compared 
to the Mjj q terms that we keep. We also are dropping terms suppressed by additional powers 
of m\. However, we keep contributions proportional to (m 2 ] /3)M 2 q logMgq, since they become 
important as we approach the chiral limit. In particular, if M 2 q < mo/3, the enhanced 
logarithm can lead to (mQ/3)logM gg ~ M qq , in which case the logarithmic term should be 
kept (it is ~ M qq ), while it is consistent to drop higher powers of (mjj/3) logM gg . In practice, 
in present simulations M 2 > ml/3, so it is likely that M qq \og(M qq ) terms are significant. 
Nevertheless, as simulations approach closer to the chiral limit, our analysis will become 
more applicable. 

A second reason for truncating the expansions is more practical. The non-analytic correc- 
tions that we have calculated are those that are determined unambiguously in terms of the 
lowest order coefficients in the chiral Lagrangian. There is no dependence of unknown scales 
~ A x , and no counterterms of the same order. This is obvious for the terms proportional 
to V and W, but is true also for those containing the X(a, b) — the scale of the logarithm 
cancels in the corrections M^ ry . This is in contrast to meson masses and decay constants, 
e.g. Eq. (|9lD, where the chiral corrections are proportional to log(M| (J /A 2 ) and thus contain 

an unknown scale. This unknown scale does enter into M^ y , but can be absorbed into the 
meson masses as described above. 

Finally, we recall that our calculation is done in the approximation that AM M qq , so 
that the octet and decuplet baryons can be treated to first approximation as degenerate. It 
would be straightforward in principle to extend our results to larger AM, although we have 
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not carried out the calculation. AM enters only in diagrams in which the internal baryon 
is an octet while the external baryon is a decuplet, or vica versa. These are the diagrams 
giving contributions proportional to C 2 . The effects of AM 7^ can be included exactly by 
shifting the mass in the baryon propagator in these diagrams, k-v — > k-v±AM. Some of the 



resulting integrals have been evaluated in Ref. [12]. For example, for the second diagram of 
Fig. (^|b), the factor of M% from the loop is replaced by an integral whose expansion for small 
AM is Ml — 3M^AMlog(M 7r )/7r. In present simulations, the second term is a relatively 
small correction. 
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